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ABSTRACT 

We consider the effect of lensing magnification on liigh. redsliift sources in tlie case tliat 
magnification varies on tlie sky, as expected in wide fields of view or within observed 
galaxy clusters. Wc give expressions for number counts, flux and flux variance as inte- 
grals over the probability distribution of the magnification. We obtain these through 
a simple mapping between averages over the observed sky and over the magnification 
probability distribution in the source plane. Our results clarify conflicting expressions 
in the literature and can be used to calculate a variety of magnification effects. We 
highlight two applications: 1. Lensing of high-z galaxies by galaxy clusters can provide 
the dominant source of scatter in SZ observations at frequencies larger than the SZ 
null. 2. The number counts of high-z galaxies with a Schechter-like luminosity function 
will be changed at high luminosities to a power law, with significant enhancement of 
the observed counts at L ^ 10 L*. 



1 INTRODUCTION 



Magnification due to gravitational lensing leads to observ- 
able effects, namely changes in the number density of galax- 
ies behind large-scale structure and galaxy clusters (known 
as magnification bias) and in the moments of the flux dis- 
tribution due to unresolved sources at high redshift. These 
and other effects of lensing magnification have been studied 
extensively in the last few decades, usually assuming simple 
expressions that apply for constant magnifications. 

In this brief note, we generalize to the case where mag- 
nification varies on the sky - the variation is taken to be 
given by a magnification probability in the source plane, 
while quantities of interested are observed as averages in 
the image plane. We apply this calculation to lensing of the 
intrinsic number count distributions of high redshift galaxies 
as well as moments of the fiux for Poisson distributed high- 
z galaxies behind galaxy clusters. Our goal is to provide 
the formulae needed for magnification effects in a variety 
of physical situations and give estimates of the scale of the 
main effects. Applications to more detailed models and re- 
sults for Sunyaev-Zel'dovich surveys have been presented in 
a separate paper (Lima, Jain & Devlin 2009). 



2 CONSTANT MAGNIFICATION 

By definition, magnification (denoted /i) is the Jacobian 
of the transformation between image (lensed) and source 
(unlensed) coordinates (e.g. Bartelmann & Schneider 2001). 
Along a given line of sight, its effect on differential solid 
angles is given by 



or fi = dQobs/dQ. We use subscript "obs" for the observed 
(or lens plane or image plane) and no subscript for the 
(unlensed) source plane. The surface brightness of galaxy 
sources, defined as the fiux per unit solid angle, is conserved 
by lensing. Since magnification increases the solid angle of 
sources by a factor n, it also increases their flux S as 



S — Sobs — l-l-S . 



(2) 



In terms of the lensing shear 7 and convergence k, the mag- 
nification is given by /i = 1/[(1 — k)^ — I7P]. 

As a result, the number density of a source population 
is modified by lensing magnification. Let dn/dS denote the 
intrinsic number density per unit fiux per unit steradian on 
the sky. Given a (constant) magnification n, it is modified 
as 



dn 
dS 



dno}^s{Sohs) 

dSohs 



1 dn 
~dS 



Sobs 



(3) 



The factor comes from transforming the angle dO, and 
the flux differential dS into their observed counterparts us- 
ing Eqns. 1 and 2. The change in argument comes from the 
fact that the observed fiux Sobs corresponds to true fiux 
S = Sobs/M- 

Given the differential number density dn/dS, we may 
define the cumulative number density n(> S), the average 
fiux of the background galaxy population per steradian S 
and the mean square fiux per steradian 5^ as 



dQ. 



dQ,o 



fidO, , 



(1) 



n{> S) = 
S = 
52 = 



—dS' 



(4) 
(5) 
(6) 
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In the presence of a constant mangification /i, the observed 
quantities are easily obtained using Eqns. 2 and 3 as: 



n(> S) ^ -n[ > 



5 



(7) 



Note that in the integrals over S for S and there is no 
upper or lower cutoff in flux. 

There is a long history in the literature of magnification 
effects on source counts (starting with Canizares 1981, 1982 
and Peacock 1982). The expressions above are consistent 
with those in the literature. We next consider the case of 
variable magnification on the sky. 



3 VARIABLE MAGNIFICATION ON THE SKY 

We wish to generalize Eqns. 3 and 7 to the case that the 
magnification varies on the sky. This variation can occur 
over a large patch of the sky with fluctuations due to large- 
scale structure or simply over the surface of a galaxy cluster 
due to variations in the surface mass density and shear over 
this surface. 

For variable magnification, the obvious step would be to 
average Eqns. 3 and 7 over the observed sky (i.e. the image 
plane), and this is indeed correct. Thus Eqn. 3 generalizes 
to driobs/dSobs = J dfi dn/dS Pobs(M)/A*^i where Pobs(M) is 
the normalized magnification probability in the image plane. 
It is often preferable to do calculations in the source plane. 
This is straightforwardly done by using the relation 



(8) 



where -P(/i) is the source plane probability. 

The above relation gives the generalized expressions: 

dnohs{Sohs)_\ _ J_ f ^^ Pjy) ^ f ■S'obs \ .g^ 

rf^obs / (m) 7 dS \ ^ J ' 

(nobB(> Sobs)) = ^ / ^M^(^) (> ^) ' (10) 



{So 

(52 



obs / — 



= s, 



(11) 

(12) 



where () denote averages of observed quantities over speci- 
fied parts of the sky. 

To obtain these results more formally, we evaluate the 
expressions on the LHS of the above equations by defining 
the average of a function X in the image plane over an ob- 
served solid angle as 



(X>ob 



(13) 



In the source (unlensed) plane the average over solid angle 
also defines P{fi) 

(X>,ourcc = ^jdQ.X = j dflPifl) X . (14) 

The function X is a function of angle 6 on the sky through 
its dependence on /i(6). Defining P(/i) on the source plane 
is conventional in lensing as it addresses questions such as, 
what is the fraction of sources that are magnified by a certain 



amount? Q Note that we have as desired 



rfM ^^P{^A = ^ M(e) = (m) = 



(15) 
(16) 



In the limit of the whole sky, we have Afiobs = Afl 
and (/i) = 1, i.e. the average magnification is unity. 

Using the relations given above in Eqns. 13-16 for an- 
gular averaging, we can obtain Eqn. 9 as follows: 



dnohs{Sohs) 

dSohs 



Af^obs 
1 



dO,o 



dllohs (Sobs 
dSohs 



(17) 



1 dn f Sobs 
d\l ^—5- -77; 



^2 dS 



— — f dfi P{fi)-— ( 

This is our first desired result. It is obviously different from 
integrating the expression for djiobs/dSobs from Eqn. 3 over 
P(/i) - doing that would have led to both factors of l//i 
being inside the integrand. Next we substitute Eqn. 3 into 



(Wobs(> Sobs)) 



1 



An„ 



dQ,ob 



dTlohs{Sohs) to/ 
777, ClJnV 



dS'o^^ 



and change variables to S' — Sobs/f* to obtain Eqn. 10. Our 
expressions for number counts agree with Schneider (2006). 
Note that the unlensed number counts are independent of 
position on the sky, as are S and S^. Eqns. 11 and 12 for the 
observed flux moments can be obtained similarly and easily 
generalize to the n-th moment as (<S"obs) ~ ifJ-"} / {l-i-) <5". 
This expression changes if there is a lower or higher limit 
to the integral over S. For instance in the case of an upper 
limit Scut, we can generalize to obtain 



<S"obs(< Scut) J ■P(^)^" ^" (< ^ 



(18) 



Applications to galaxy clusters are discussed below. An- 
other application is the contribution of unresolved point 
sources to CMB anisotropies, given by d — S^{< Scut). 
The upper cutoff Scut is usually introduced to remove re- 
solved objects brighter than the cutoff. With lensing the 
observed contribution is enhanced - given by the above re- 
sult with n = 2. The enhancement depends on the slope of 
the dn/dS relation at the cutoff. Finally we note that with a 
flux limit Eqn. 11 is no longer true, since surface brightness 
is only conserved when integrated over all fluxes. 



4 GALAXY CLUSTERS 

Galaxy clusters produce magnifications ranging from ~ 10% 
enhancements above unity to factors of several or more as 



^ One must of course ensure that theoretical predictions are also 
in the source plane. This occurs naturally in ray tracing simula- 
tions which start the rays at the observer and trace backwards 
(e.g. Jain, Seljak & White 2000). However predictions that rely on 
the Born approximation apply to the image plane. See Hilbert et 
al. (2008) for a discussion of simulation predictions; they also con- 
sider the effects of multiple imaging which we are not concerned 
with here. 
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Figure 1. Left: Mean magnification (fj,) and mean squared (^^) measured within clusters at redshift 2 = 0.5 as a function of cluster mass. 
The upper panel averages within the virial radius of the cluster while the lower panel uses a smaller radius within which the overdensity 
is 2500 times the critical density. In X-ray and SZ analysis, different choices for the cluster radius are made which typically lie within 
these two. Right: Contamination of submillimeter galaxies to the SZ flux within the virial radius of clusters of virial mass M^ir = lO^'*'^ 
and 10^^/i~^Mq. The solid line shows the intrinsic SZ flux AS^^. The light shaded region delineate the 1-a region from contamination 
of submillimeter galaxies from Poisson noise in their counts, i.e. tr^^j = S^vir- The dark shaded region accounts for lensing magnification, 
which further enhances the noise by (fj,^) / {fj.) , cf. Eqn 22. 



one approaches the critical curves. As a resuh both num- 
ber counts of background galaxies and the flux moments of 
unresolved background sources are significantly altered. 

Consider the unlensed average number of background 
galaxies within a cluster solid angle AJlvir defined by its 
virial radius, i.e. A^'vir = Af2vir n{> S). Notice that in the 
unlensed case Af2vir ~ Af2, i.e. the virial radius of the clus- 
ter is actually the intrinsic angle. With lensing the virial 
radius is now the observed solid angle, i.e. Af2vir ~ ASlobs, 
and we have 



(iVobs(> S'obs)>vir = Af^vir(nobs(> Sobs)} (19) 



Note that is now the magnification probability within 

the cluster virial radius (in the source plane as before). In 
the limit of constant magnification inside the virial radius 
we get Aobs = N/fi. Furthermore, if there is no lower limit 
to the flux integral Sobs — >■ 0, we have {Nobs) = N/{fi), i.e. 
we observe fewer galaxies in the line of sight of clusters. Note 
that our result for Nobs appears to differ from some of the 
literature (e.g. Schneider 2006), but the difference is that we 
use the same solid angle for the observed and unlensed case, 
because the only solid angle in town is the observed size of 
the galaxy cluster. 

The intrinsic mean flux and mean square flux within 
the cluster solid angle Ailvir are 5vir = Afiyir S and = 



Ailvir S^. With lensing we have 

(*5obs)vir = Ai^vir{*5obs) = *5vir , (^0) 
(*5^obs)vir = Ar^vir (5^obs) = ~^ — y^'^vir- (21) 

Consider a cluster of observed radius 1 arcmin. The 
average flux from background galaxies crossing this 1 arcmin 
cluster is the same with or without lensing: with lensing, the 
background galax;ies are brighter but less numerous and the 
two effects exactly compensate. So, lensing does not create a 
bias by increasing the expected flux within the cluster solid 
angle relative to a random 1 arcmin patch of the sky, as long 
as there is no flux cutoff. 

For SZ surveys, the mean flux from background galaxies 
is subtracted in obtaining the SZ decrement or increment 
(depending on the observed frequency), but there is addi- 
tional scatter due to the Poisson error from shot noise fluc- 
tuations (a given cluster having more or less galaxies than 
the expected average). This effect is important for SZ sur- 
veys, where high-z submillimeter galaxies can contaminate 
the SZ signal (White & Majumdar 2004; Knox et al. 2004; 
Lima, Jain & Devlin 2009). Lensing enhances this source of 
scatter signiflcantly, as the factor (/x^ for galaxy clus- 
ters can be quite large - see Fig. [T] The estimates shown are 
based on an analytical model of cluster halos that uses NFW 
profiles and elliptical iso-potential contours (Lima, Jain & 
Devlin 2009). Our estimates are conservative in that we in- 
clude cluster ellipticity but not substructure. Fig. 2 also in- 
dicates the error made in approximating moments of the 
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Figure 2. Left panel: Intrinsic and lensed galaxy number density dn/dlnS using different lensing prescriptions. The true underlying 
distribution is assumed to be of a Schechter type, which is then lensed by intervening halos. The solid (green) curves use Eqn. 9; the 
two alternatives shown (and described in Section 5.2) can change the predicted counts by an order of magnitude or more at the high 
luminosity end. Right panel: As above, but for different source redshifts. 



flux by using the mean magnification (e.g. Refregier & Loeb 
1997). 

Notice that if we attempt to remove galaxies above some 
flux Scut, which could be removed if they were resolved, we 
introduce a difference between the observed and intrinsic av- 
erage flux through the cluster, i.e. (5obs) 7^ S. In that case, 
the mean CMB flux is not equally contaminated by back- 
ground galaxies and subtracting this mean flux from the 
cluster flux does not cancel the galaxy contribution on av- 
erage. Therefore, removing bright galaxies from the sample 
biases the SZ signal. An observationally relevant situation 
arises if a flux cutoff is used to identify sources - thus al- 
tering both the number counts and flux contribution from 
unresolved sources (e.g. Refregier & Loeb 1997). 

5 ILLUSTRATIVE NUMBER-MAGNITUDE 
RELATIONS 

5.1 Power law 

It is common in the magnification bias literature to consider 
the case of a local power law in the logarithm of the number 
counts. We then obtain for the observed cumulative number 
density: 

n(> S) (X nob3(> Sobs) = ^ n{> Sobs). (22) 

W 

Working with apparent magnitudes instead of fluxes this 
gives (using m = — 2.51og]^gS -I- Constant) 

n{< m) (xm" — )> nobs(< mobs) = , , n(< mobs). (23) 

Note that the use of the power law for n{> S) allows us to 
simplify the integral over /x. The above equation agrees with 



the standard expression (e.g. Broadhurst, Taylor & Peacock 
1995): nobs(< m) — fi '^""^ n(< m) for the case of con- 
stant magnification. For variable magnification, one needs 
to evaluate the averages as above. If one simply uses the 
factor {fi'^'^"~^) behind a galaxy cluster instead of Eqn. 1231 
an error in the number counts can result. The error ranges 
from 3% for a cluster of mass lO^'' /i~^Mq to 6% for a mass 
of 10^^ h~'^MQ. This would result in an equivalent bias in 
the inferred cluster mass. Mass estimates that rely on num- 
ber counts may be feasible for large samples of clusters from 
future surveys. 

To contrast the results from Eqn. 9 with other formulae 
used in the literature, consider first the naive generaliza- 
tion to variable magnification (which is correct in the image 
plane) : 

„ .... dn-obs(Sobs) f , -P(m) dn / Sobs \ ,„.^ 

Prescription A : — ^ '- = / d^^^f^ — . 24 

rfSobs J y? dS \ ^. J 

This underestimates the lensing effect. Alternatively some 
authors drop the /i factors altogether (Paciga, Scott & 
Chapin 2009), which overestimates the lensing effect: 

5.2 Schechter luminosity function 

A single power law number-magnitude relation is lensed into 
an observed relation with same power law (but different am- 
plitude). However if the intrinsic dn/dS is not a power law, 
then lensing changes the shape of the distribution as well. 
High magnification events shift galaxies with low fluxes to 
high fluxes - hence if dn/dS falls sharply at high S, magni- 
fication can significantly enhance the counts at these fluxes. 



© 0000 RAS, MNRAS 000, 000-000 



Magnification Effects on Source Counts and Fluxes 5 



We illustrate the effect of magnification for realis- 
tic galaxy populations by considering a Schechter function 
(Schechter 1976): 



dn{S) 
dS 



-s/s' 



(26) 



For galaxies observed in a narrow redshift interval, such a 
dn/dS relation can arise due to the Schechter luminosity 
function of the population. 

In Fig. [2] we show the intrinsic distribution and the 
lensed versions according to the two incorrect prescriptions 
(denoted A and B) mentioned above, as well as the cor- 
rect prescription of Eqn. 9. We assume all galaxies are at 
redshift Zs — 2.1 and use a P(pi) obtained from N-body 
simulations by Hilbert et al (2007). We also choose a — 
in the Schechter function. We plot dn/d\nS to relate more 
easily to the observational literature which shows number 
per unit absolute magnitude. Our results in Fig. 2 may be 
matched with high-z luminosity functions by replacing S/S* 
with L/L*. 

The lensing contribution (correctly included in the 
green curves) is large for the Schechter function at the bright 
end {S/S* > 10). Thus for a population of galaxies with a 
Schechter luminosity function, the observed dn/dS will not 
retain the exponential tail of the Schechter function. Using 
Eqn. 9 and approximating the Schechter function as having 
a sharp cutoff at 5 = 5**, it is easy to see that the integra- 
tion over generates a power law in dn/d\nS whose slope 
is —2, due to the asymptotic slope P(/i) oc The green 
solid curve approaches this slope beyond S — lOS* . 

Choosing a lower (higher) value of a slightly enhances 
(suppresses) the lensing contribution at fixed S/S*. Fig. [2] 
also shows that the difference between the three prescrip- 
tions is large at the bright end for the Schechter function. 
The dn/dS in Fig. 2 is normalized so that the distribution 
matches submillimeter galaxies measured by the Balloon- 
borne Large Aperture Telescope (BLAST) (Devlin et al 

2009) at 500 /im. In a separate study we show that if the 
sub-mm galaxies lie at z ~ 2 — 3, then lensing of an in- 
trinsic Schechter function can explain the observations from 
BLAST and other surveys (Lima, Jain, Devlin & Aguirre 

2010) . 

Similar comparisons to luminosity function observations 
in the visible bands can be carried out: it is especially im- 
portant to include lensing in the comparison of \ow-z data 
with data at z ^ 1 since the magnification contribution is 
significant only at high-z. 



The quantitative estimates presented here for galaxy 
clusters are based on analytical models of halos. While these 
incorporate realistic density profiles as well as halo elliptic- 
ity, they miss the full complexity of halo bimodality (due to 
major mergers) and substructure. These features only en- 
hance the effects of averaging we have considered. 

Finally we note that the lensing effects discussed here 
do not impact the predicted (magnification induced) cross- 
correlation of number counts measured in different redshift 
bins (Moessner & Jain 1998) . Such cross-correlations depend 
on the two-point cross-correlations of magnification with the 
galaxy density. Thus measurements of magnification bias 
from galaxy-quasar cross-correlations or of its contaminat- 
ing effect on high-z ISW cross-correlations are unaffected by 
the spatial averaging issues discussed here. 
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6 DISCUSSION 

We have derived expressions for computing average quanti- 
ties in the observed image plane, given a distribution P(/i) 
of magnifications in the source plane. Our results, summa- 
rized in Eqns. 9-12, generalize expressions found in the lens- 
ing literature for the case of constant magnifications. We 
illustrated the effect of lensing on steep number counts of 
background galaxies and on boosting the contamination that 
high redshift galaxies induce in cluster SZ fluxes. The for- 
mulae we have presented may be useful for studying the 
intrinsic properties of high redshift galaxies and for current 
and upcoming cluster surveys. 
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